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ABSTRACT 



Laminar natural convection at unity Schmidt number over a horizontal surface with 
a weak normal velocity at the wall is studied using an integral analysis. To charac- 
terise the strength of the blowing, we define a non-dimensional parameter called the 
blowing parameter. After benchmarking with the no blowing case, the effect of the 
blowing parameter on boundary layer thickness, velocity and concentration profiles 
is obtained. Weak blowing is seen to increase the wall shear stress. For blowing 
parameters greater than unity, the diffusional fiux at the wall becomes negligible 
and the flux is almost entirely due to the blowing. 



Introduction 



Fluid motion driven by an indirectly induced pressure gradient normal to the direction of the 
density potential difference is termed as indirect natural convection P; 0] . An example is natural 
convection occurring over a horizontal heated surface. At large Grashoff numbers, the flow is of 
boundary layer type , and can be described by similarity solutions Such boundary layers 

could form over a porous surface with blowing at the surface, the effect of which on the behavior of 
the boundary layer is not obvious. Previous investigations on the effect of blowing on indirect free 
convection boundary layers were conducted by Clarke & Riley j^] and Gill et.al. [j] using similarity 



1 



analysis of the boundary layer equations. The boundary layer equations admit similarity solutions 
only for a restrictive and interdependent spatial distribution of wall velocity and density potential. 
Hence, the scope of these solutions is limited. We study indirect natural convection boundary 
layers for unity Schmidt(or Prandtl) number fluids, when a weak normal velocity is imposed at 
the wall. The wall velocity and density potential are not functions of horizontal distance. The 
imposed velocities are considered to be weak in the sense that they are of the order of vertical 
characteristic velocity in the boundary layer when no blowing is present. An integral method is 
used to obtain approximate solutions to the boundary layer equations for blowing velocities and 
wall concentrations that are constant in the horizontal direction. 



The notation used is plain symbols for nondimensional variables, superscript ~ for dimensional 
variables and subscript c for characteristic scales. The paper is organised in the following manner. 
Using the characteristic scales in the problem, the boundary layer equations are derived. Integral 
equations are formulated by integrating the boundary layer equations across the boundary layer 
thickness, using profiles that satisfy the boundary conditions. The solution procedure of the re- 
sulting initial value problem is exp ained and the results are initially compared with the similarity 
solutions of Rotem and Classen [4| for no blowing case. The effect of blowing on the boundary 
layer thickness, concentration and velocity profiles, and mass flux are studied. 



Boundary layer approximation 
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FIG. 1 

Schematic of the problem studied 



Consider a fluid having a concentration Co of some species above a horizontal porous surface 
maintained at C^, with Cw < Co- Also, a fluid with concentration C^j flows through the surface 
with a normal velocity Vw (The problem is equivalent to natural convection over a heated plate 
with blowing, with Cw representing the wall temperature and Co the temperature of the ambient.) 
The geometry of the problem studied along with the variables used are shown in Figure u, v 
are the velocity components in the x,y directions, and the boundary layer thickness is (5(x).The 
wall velocity is considered to be independent of the density potential. 

Let the characteristic length scales in the horizontal and vertical directions be Xc = L, and 
He = ^c, a characteristic boundary layer thickness. L is taken as the horizontal dimension of the 



surface. Let Ucandvc denote the characteristic horizontal and vertical velocities, Pc the character- 
istic pressure and Cc = Cq — Cw = AC* the characteristic concentration difference. 



Following Rotem and Classen for 5*0 ~ 1 

~ , «c ~ jGrT , ~ ^Gr]!\ ~ p{u/LfGrf (1) 

Gr^ ^ ^ 



where, Gr^ = ^^^^ , the Grashoff number. Here g,(3,v denote acceleration due to gravity, 
salinity volumetric expansion coefficient and kinematic viscosity respectively. The horizontal char- 
acteristic velocity scale can be rewritten as. 



g(3AC6c , Vc^Uc {Sc/L) , pc ~ pul (2) 



Uc can hence be viewed as a free fall velocity over the characteristic boundary layer thickness. 
The vertical characteristic velocity is 6c/L times smaller than Uc and the predominant balance is 
between motion pressure and the horizontal velocities. The motion pressure, p = ps — Ph, where 
Ps is the local static pressure in the presence of density variations, and p^ is the local hydrostatic 
pressure when the fluid through out is at the reference density p^- 

Normalising the governing equations of two dimensional continuity, Navier-Stokes and species 
equation with these characteristic scales, dropping terms of order < l/Grf^^ when Gri — > oo, and 
assuming Schmidt number, ~ 1 we obtain the non-dimensional boundary layer equations as, 

du dv ^ 
dx dy 



du du dp d'^u 

dx dy dx dy"^ 

dp 
dy 



^-.+-- = --. + 1^ (4) 



(5) 



dc dc d'^c , , 

dx dy dy'^ 



with the boundary conditions. 



du dc 
at y =5{x) : n = 0, — = 0, c = 0, — = (7) 

aty = 0: u = 0, c=l v = —= = K., (8) 

Vc GrY 

where x = x/L, y = yjbc, u = u/uc, v = vjvc, p = p/pc, c = 6{x) = 5{x)/6c, Rcu, = 

VyjL/u. Equation shows that the vertical concentration distribution creates a vertical pressure 
distribution, the horizontal gradient of which in Q drives the flow. The concentration distribution 
is, in turn, a result of this motion as well as diffusion. 



The Blowing parameter is the non-dimensional normal velocity at the wall. The blowing 
parameter also shows the condition for similarity of the boundary layers in the presence of a weak 
blowing. When the normal velocity at the wall has a dependence of ~ x"^/^ the boundary 
condition become independent of x and a similarity solution can be obtained. The normalised 
terms in the boundary layer equations are of order one. Hence these approximate equations can 
be expected to hold till the blowing parameter -^^g- ~ 0(1). We restrict our calculations up to 

= 3. 



Integral Formulation 



Consider a steady laminar indirect natural convection boundary layer with a weak blowing that 
could be described by the boundary layer equations © - (©• Integrating ^ from the surface to 
5{x) and applying the boundary conditions, we get the non dimensional integral mass conservation 
equation as 



dx 







n dy + t>5 - K, = (9) 



where, vs = v^/vc is the non-dimensional entrainment velocity at the edge of the boundary layer. 
The non-dimensional integral x momentum equation is obtained by integrating the equation ob- 
tained after multiplying (jH)) by u and subtracting Q). 

d r^^^') f^(^) r)r) fill 

^ ^'dy+ ^dy + p =0 (10) 

dx Jo Jo dx dy\^^^ 

This equation expresses a balance of horizontal momentum flux with horizontal pressure gradient 
and wall shear stress. Integrating Q from to 5{x) and substituting for vs from Q, the non- 
dimensional integral species equation is obtained as 

d f^^'-^"^ dc 
_« / ucdy + ^ +V^ = (11) 
ax Jo '^y\y^Q 

The terms in the above equation represent convection by the horizontal velocity and entrainment 
from ambient, the diffusive flux at the wall and the blowing flux at the wall respectively. The 
concentration profile is approximated by the function, 

/(r?,a) = c =(l-r?)2 (1 + ar?) (12) 

where, r/ = y/6{x) and the coefficient a captures the non similar nature of the profile in the 
presence of blowing. The profile H12() satisfies the boundary conditions Q and © ie., /(0,a) = 
1, /(I, a) = 0, /(I, a) = 0, where, superscript . denote differentiation with respect to rj. 



The coefficient a in (fT^ is found from the condition that the profile satisfies the species equation 



^ at the wall ie, 

/(0,a) =K,5/(0,a) (13) 
Substituting the values of /(O, a) = a — 2 and /(O, a) = 2 — 4a in (fT^j) and solving for a,we get 
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a(F^,x) = 2- 

4 + VV; 

The velocity profile is approximated by the polynomial 

u{r], b) = ci + b{x) T] + C2rf' + C'irf + C4 rf 



(14) 



(15) 



where the dimensionless coefficient h captures the non similar nature of the profile in the presence 
of blowing. We need to satisfy the boundary conditions (0) and 0, ie. n(0, 6) = 0, m(1,6) = 
0, andn(l,6) = 0, and the condition obtained from the x-momentum equation at the wall, 



— u(0) - (dh) - = 



(16) 



where ' denotes differentiation with respect to x. In the above equation, the horizontal pressure 
gradient in Q) is expressed in terms of the concentration profile using © as 



dp 
dx 



r)=0 



-{dh)' where, Ii = / f{vi)dri 
Jo 



(17) 



Using the above conditions, ci = 0, and C2, C3, C4 can be expressed as functions of 6, Vw, and 5. 
The final expressions of the coefficients are given in the Appendix. The velocity profile is obtained 
as 



u{r],b) ^ ri{ri- if 



(2 + r?(4 + y^(5)) - 



(18) 



2 ' ■ 4(4 + K, 5)2 

Rewriting (lllj) and (|1())) in terms of the profiles (|12|) and (|15)) , we get the integral equations for 
species and momentum. 



[6 lufl -Vw-\ — = 



[542 -6^h]' + 



u{0,b) 



where. 



luf = / u{7], x) f{r], a) di], /„2 / u{r], bf di]; L 



/(ry, a) di] di] 



(19) 
(20) 

(21) 



In (|2flj) above, the horizontal pressure gradient in (|10|) has been replaced in terms of the concen- 
tration profile by integrating the y momentum equation ^ from 1 to as below 

r^(^) Pin 

Substituting the values 42, 4, n(0), /(O) and simplifying, we get two autonomous nonlinear 
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FIG. 2 

Comparison with Rotem and Classen 0| for no blowing case, (a) Non dimensional boundary layer 
thickness (b) Non-dimensional concentration profile (c) Non-dimensional Velocity profile 

second order ordinary differential equations for the unknowns 5{x) and b{x) as follows. 

5" + 5' {A6' + B) + Cb' + D = (23) 
6" (1 + E6') + 5' (^Pf + G6' + + h' {15' + J) + K = {) (24) 

where the coefficients, A to K are nonlinear functions of 6{x), Vyj and6(x) The expressions for these 
coefficients are given in the Appendix. 

Equations H23|) and (|24() needs three initial conditions (5(0), (5'(0), and 6(0) for solution as 
an initial value problem. We assume that as x — > the boundary layer behaves similar to the 
case with no blowing. This is equivalent to assuming an infinitesimal distance where no blowing 
effects are present. Hence, we use the values of 5(e) = 0.0179, 5'{e) = 6369.7 and 6(e) = 0.2358, 
where e = 10~^, as calculated from the similarity solutions of Rotem and Classen |J|as the initial 
conditions. It was found that the solution was insensitive to changes in initial conditions up to an 
order magnitude from the above conditions. The NDSolve routine of MATHEMATICA was used 
to numerically solve (|23|) and (|24|) simultaneously. 



Results and Discussion 



The results are initially compared with similarity solutions for the no blowing case by Rotem 
and Classen 



|. The integral analysis gives 5{x) ~ x^/^, same as that obtained in the similarity 
solution. The boundary layer thickness 6 in the integral analysis is equal to the thickness at which 
c = 0.05 in the solution of Rotem and Classen Q]. Figure (2(a)) compares the x dependence of 
the non dimensional boundary layer thickness 5 with that of Rotem and Classen Q]. Figures (2(b) ) 
and ( 2(c) ) compare the concentration and velocity profiles obtained from the integral analysis with 
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FIG. 3 

Effect of blowing on (a) boundary layer thickness and (b) concentration profile 

those from the similarity solution. Considering the approximate nature of integral methods, the 
comparison is satisfactory. 



Figure (3(a) ) shows the effect of blowing on the boundary layer thickness. As expected, higher 
blowing results in larger boundary layer thickness; the increase in total added fluid increases the 
boundary layer thickness. The concentration profiles at x = 0.5 for different blowing parameters 
are shown in Figure ( |3(b)| ). The reduction in the gradient of concentration at the wall shows that 
the diffusive flux decreases with increase in the blowing parameter. 

The major effect of the change in concentration profile shape is on the pressure distribution. 
The motion pressure outside the boundary layer is zero. The motion pressure within the boundary 
layer is given by p{y,x) = —pgfi j^~^{Co — c)dy. Presence of increased amounts lighter fluid 
in the boundary layer due to blowing results in larger value of Co — c. This results in a larger 
favourable horizontal pressure difference with the outside motion pressure than in the case of 



without blowing, resulting in larger velocities. This could be observed in Figure ( 4(a) I shows the non 
dimensional horizontal velocity distribution for different blowing parameters at x=0.5. Figure p(b)| 
shows the variation of horizontal pressure gradient at the wall as a function of horizontal distance for 
increasing value of the blowing parameter. The pressure gradient was calculated using the relation 
^ = —pgP-^ iif=o ^i^o — c)dy. The increased amounts of lighter fluid in the boundary layer has 
the effect of making the horizontal pressure gradient at the wall more favourable. This increases ^ 
at the wall and results in increased wall shear stresses with blowing. This behaviour is in contrast 
to what is observed in a shear boundary layer with blowing. For example, in the case of the Blasius 
boundary layer, blowing reduces the wall shear stress. The vertical characteristic velocity in the 
boundary layer is ^ times the horizontal characteristic velocity, as can be seen from Equation 0. 
Therefore, the vertical momentum from blowing velocities of the order oivcie, ~ 0{1), doesn't 
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FIG. 4 

Effect of blowing on (a)liorizontal velocity profile and (b)horizontal pressure gradient at the wall 

seem to be sufficient to cause a change in the near wall linearity of the horizontal velocity profile. 
For the current range of blowing parameters investigated, the major effect seems to be due to 
increase in motion pressure, which increases the horizontal velocities, and the wall shear stress. 



The total mass flux of the species into the boundary layer is 



dc 

VrnAC + D—' 



y=o 



(25) 



The first term on the right hand side of ()25() is due to blowing(advection) and the second term 
represents diffusion. The non-dimensional flux is. 



Or 



(26) 



The variation of non-dimensional diffusive flux over the length of the boundary layer is shown 
in Figure (5) b Blowing reduces the diffusive flux. Note that when Vw > 1, except very near the 
leading edge, the diffusive flux is an order lower than the blowing flux, which is equal to V^. 

Figure © shows the longitudinally averaged blowing and diffusive fluxes as a function of Vw 
The total flux is also shown in the flgure. As pointed above, the mean diffusive flux drops with 
Vuj- At Vw = 3, the diffusive flux is about ^rd the diffusive flux at no blowing. The blowing flux 
increases linearly with the blowing parameter, and at Vw = 3, the blowing flux is about flfteen 
times the diffusive flux, and about 5 times the no blowing flux. The variation of the ratio of the 

total flux with blowing to the flux without blowing {-frfoy) could also be noticed from Figure ©. 

s 

For example, the flux with Vu,=l is double the flux with no blowing. To take a speciflc example. 
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Variation of the longitudinal distribution of Variation of non-dimensional fluxes with blowing 
diffusive flux with blowing parameter parameter 



this means that at a commonly encountered Gr^ = 10^ in water, blowing velocities of the order 
of 1 mm/s doubles the flux from the no blowing case. Hence, weak blowing increases the flux 
appreciably over the corresponding no blowing flux. 



Conclusions 



The integral analysis developed in this paper is used to study the effects of weak blowing on 
indirect laminar natural convection boundary layers. The blowing velocity and wall concentration 
are constant along the wall. The analysis is subject to the assumptions of small (5/-L (large Grashoff 
numbers), unity Schmidt(or Prandtl) numbers and blowing velocities of the same order as the ver- 
tical characteristic velocities due to free convection in the boundary layer. Approximate solutions 
were obtained by numerically solving the integrated boundary layer equations. Following are the 
major conclusions from the study: (a) As expected, the boundary layer thickness increases with 
blowing (Figure 3(a) ). (b) Weak blowing increases the horizontal velocities in the boundary layer, 
resulting in the wall shear stress also increasing with blowing (Figure 4(a) ). (c) The concentration 
gradient at the wall and thus the diffusive flux reduces with increase in the blowing parameter 
(Figure [3(b)| ). (d) For Vw > 1, the flux due to advection caused by blowing becomes the predomi- 
nant contribution to the total flux and hence the total flux shows a linear dependence on Vw (e) 
Comparison of total flux with that in the case of no blowing show that a weak blowing results in 
considerable increase in flux (Figure inj. 



Appendix 



The coefficients in the velocity profile p5j) are, 

c2 = 56 b- Qd^6', c3 = -b (3 + 106) +2QS^ S', c4 = ^^^^^iMl _ g ^2 y ^27) 

(2 + 10 5) (6 + 10 5) 

where, Q = — ^ — - (28 

4(4 + 105)2 

The coefficients in the differential equations (|23|1 and (|24j) are obtained as 



3/1 3 1 7 ^ 

^ = ^ + (^2 + K,5 ~ 4 + V^5 ^ 6 + Vy,d ^ W + TK^J ^^^^ 
-2b i4 + V^6) (816 + 7K,5 (104 + K, J (25 + 2^;^))) 

53 (2 + K, 5) (6 + K, S) (19 + 7K; ^ ^ 

-2(4 + K,(5)' (204 + 7K;5 (13 + K,5)) ^ 1680 (4 + 5f (6 + ¥^5(4 + 6)) 
5^ (2 + <5) (6 + 6) (19 + 7V^6) ' 5^^ [2 + 5) (6 + 6) (19 + 7K« 5) ^ ^ 



,5^ (2 + S) (6 + 5) p_ -.5 (240 + J (236 + 5 ,5 (12 + K, ,5))) 
'bii + V^5f {17 + 2V^5y ~ 2b(^4 + V^5f {17 + 2V^5) 



G=^ + V (—^ 2 1 2 \ 

"'V2 + K;5 4 + K;(5 6 + K;J 17 + 2145/ ^ ^ 



504(5 (2 + K;(5) (16 + 3K;<5) 
(2 + <5) (6 + -5) (17 + 2K, ,5) 
2 6^ (4 + K,(5)2 (76 + -14 ,5 (34 + 3^^(5)) 
■ 6^3 (2 + K,5) (6 + 14(5) (17 + 2K,(5) 



(34) 



J 



-4 (4 + '^)' (76 + K; (5 (17 + 14 ^)) 
(52 (2 + 14 S) (6 + 14 S) (17 + 2 14 (5) 

-5040(4 + 14(5)2 
^ (54 (2 + K;,5) (6 + 14 <5) (17 + 214 <5) 



(35) 
(36) 
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